Pomeranchuk cooling of the SU(2A^) ultra-cold fermions in optical lattices 
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We investigate the thermodynamic properties of a half-filled S\J{2N) Fermi-Hubbard model in 
the two-dimensional square lattice using the determinantal quantum Monte Carlo simulation, which 
is free of the fermion "sign problem". The large number of hyperfine-spin components enhances 
spin fluctuations, which facilitates the Pomeranchuk cooling to temperatures comparable to the 
superexchange energy scale at the case of SU(6). Various quantities including entropy, charge 
fluctuation, and spin correlations have been calculated. 
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The SU(2A^) and Sp{2N) Fermi-Hubbard and Heisen- 
berg models were originally introduced to condensed 
matter physics as a mathematic convenience 1-4]. They 
generalize the SU(2) models of electrons and enable the 
1/iV-expansion to systematically handle strong correla- 
tions. In the cold atom context, it was proved that the 
large spin Fermi systems with the hyperfine spin F — ^ 
possess a generic Sp(4), or, isomorphically S0(5) sym- 
metry, without fine tuning Q . It can be augmented 
to the even higher symmetry of SU(4) with alkaline-earth 
fermions due to their fully filled electron shells f3] . In this 
case, their hyperfine spins completely come from nuclear 
spins, thus interactions are spin-independent. The recent 
development of ultra-cold Fermi gases has provided an 
opportunity to realize these high symmetries, thus they 
are not just of academic interests. Experimentally, the 
SU(6) symmetric system of ^^■^Yb with F = I ^ ^ and 
the SU(IO) symmetric one of ^^Sr with F = / = | have 
been cooled down to quantum degeneracy [T-Qj. The 
study of SU(2iV) by using alklaine-earth fermions was 
also proposed in Ref. (lo| . 

The properties of large spin ultra-cold Fermi gases with 
SU(2iV) and Sp(2iV) symmetries are fundamentally dif- 
ferent from those of large spin solid state systems . In 
the former case, quantum spin fluctuations are enhanced 
by the large number of fermion components denoted by 
27V — 2_F-|- 1, which are even stronger than those in spin- 
i systems. On the contrary, the large value of spin sup- 
presses quantum fluctuations in solid state systems, and 
drives the system to the classic limit. The novel physics 
in large-spin ultra-cold Fermi systems with high symme- 
tries includes unconventional superfluidity |12h16I |. ex- 
otic quantum magnetism [17H2ll |. In particular, various 
SU(27V) valence-bond solid and spin liquids have been 
proposed that have not been observed in solid state sys- 
tems [ii-il. 



In current experiments, ultracold fermions in optical 
lattices have been cooled down to temperatures below the 
hopping energy scale of t, but, above the spin superex- 



change energy scale of J yet. This motivates us to study 
the thermodynamic properties of the ultracold fermions 
in this temperature regime. On the other hand, the 
strong correlation physics of the Fermi-Hubbard model 
in this intermediate temperature regime is also interest- 
ing, which has not been fully explored yet in solid state 
physics context. In cuprates, J is around lOOOK at which 
lattices melt. In this temperature regime, the intertwined 
effect between the quantum and thermodynamic fluctu- 
ations leads to interesting phenomena. The quantum 
critical phenomena have been studied in the SU(2) case 
[m 01- Recently, the high temperature properties of 
SU(2A^) Hubbard model has been studied from series ex- 
pansions, which is only accurate at T ^ max(t, U) [27|. 

In this paper, we focus on the intermediate temper- 
ature regime {t > T J) of current experiment in- 
terests and study the thermodynamic properties of the 
half-filled SU(2A^) Fermi-Hubbard model in 2D square 
lattice using the determinantal Quantum Monte Carlo 
(DQMC) simulation [H,!!^, which is an unbiased, non- 
perturbative method. It is free of the sign problem at 
half-filling, which enable us to obtain reliable results 
down to low temperatures {T/t ~ 0.1). Various ther- 
modynamic quantities are calculated for different values 
of 2N. Special attentions are devoted to the interaction- 
induced adiabatic cooling, that the system is cooled down 
to the temperature scale around J by adiabatically in- 
creasing interactions. This Pomeranchuk cooling effect, 
though it is very weak in the SU(2) case 3^, 3l|, appears 
in the SU(6) case, which can cool the system to temper- 
ature lower than J. This is relevant to the experiment 
systems of ^^'^Yb. 

We consider the following SU(2A^) Fermi-Hubbard 
model defined in the 2D square lattice at half-filling as 

U 
2" 

where a runs over the 2N components; denotes 
the summation over the nearest neighbors; Ui is the to- 



H--tJ2 + h.c.) + -J2{n,~ nY, (1) 
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tal particle number operator on site i defined as = 
X]a=i '^L'^ia- Eq.(IT|) is invariant under the particle-hole 
transformation in bipartite lattices. Similarly to the case 
of SU(2), it is easy to prove that the sign problem is also 
absent for the half-filled S\J(2N) Hubbard model of Eq. 
[T]in bipartite lattices in the DQMC simulations. 

Below we will present our DQMC simulations for ther- 
modynamic quantities for the SU(2iV) Hubbard model 
with 27V = 4 and 6 in the square lattice oi Lx L with the 
periodical boundary condition. The second order Suzuki- 
Trotter decomposition is used. The Trotter steps are 
taken At — f3/M where (3 = l/T is the inverse of the 
temperature T and M ranges from 30 to 100 depending 
on temperatures. We have checked that the convergence 
of the simulation results with respect to different values 
of At. Instead of using the Hubbard-Stratonovich (HS) 
transformation in the spin channel [32,] , we adopt the HS 
transformation in the charge channel which maintains the 
SU{2N) symmetry explicitly [s^. This method gives rise 
to errors of the order of (Ar)^. 

We will first address the possibility of cooling down the 
system by adiabatically increasing the repulsive interac- 
tion. For spinful fermion systems (e.g. '^He), the Pomer- 
anchuk effect refers to the fact that increasing tempera- 
tures can lead to solidification because the entropy (per 
particle) in the localized solid phase is larger than that 
of the itinerant Fermi liquid phase. The reason is that 
in the Fermi liquid phase, only fermions close to Fermi 
surface within T contribute to entropy, while in solids 
each site contributes to nearly ln2 « 0.69 if T is com- 
parable to the spin exchange energy scale of J, which is 
much smaller than the Fermi energy. J/t « At/U in the 

strong coupling regime, and J/t cx e^'^'^^f^ in the weak 
coupling regime. It reaches maximum in the intermedi- 
ate coupling regime. In the lattice systems near or at 
half-filling, increasing interaction suppresses charge fiuc- 
tuations and drives systems to the Mott-insulating state, 
thus we would expect Pomeranchuk cooling during adi- 
abatically increasing interactions [33|. However, for the 
SU(2) Hubbard model, both at 2D and 3D, DQMC sim- 
ulations show that the effect of Pomeranchuk cooling is 
not obvious at half-filling [s^, Isll, [sl] . This is because for 
the SU(2) case, the tendency towards magnetic ordering 
is strong at temperatures lower than J, which suppresses 
the entropy capability in the Mott insulating state (30| . 

In Fig. [Jl curves of the simulated entropy per parti- 
cle (not per site) S^^iiN) 2A^, are plotted. Ssu(2N) 
is defined as Ssu(2N) = S/{NL'^), where S is the total 
entropy in the lattice. It is calculated from the formula 



,{2N) 



(T) 



ln4 



E{T) 
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where ln4 w 1.38 is the entropy at the infinite tempera- 
ture, or, equivalently, T ^ U; E{T) denotes the average 
internal energy per particle at temperature T. The inter- 
action parameters are chosen for the free case {U /t = 0), 
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FIG. 1: Entropy per particle Ssu(2N) for the SU(2A'') Hubbard 
model at half-filling v.s. 2A'' in a 10 x 10 square lattice. The 
temperature is fixed at T/t — |. The line of C//t = oo is from 
the results of Eq. (3] 



the intermediate {U/t = 8) and strong {U/t — 12) inter- 
action strengths, and also the limit oi U/t = oo. The 
temperature T/t = is low compared to the band width 
W/t = 8, and is comparable to J for U/t = 8 and 12 
around which spin exchange is strongest. The results of 
the SU(2) case are consistent with those in Ref. [sof . 
At U/t = 8 and 12, Ssu{2) becomes even smaller than 
that at U/t — 0, which shows effectiveness of suppressing 
entropy from spin correlations. Consequentially, we do 
not expect the Pomeranchuk cooling in this interaction 
regime. If further increasing U/t to very large values, 
finally J will become much smaller than T, and in this 
case we expect that Sgu{2)/kB will grow and saturate at 
the In 2 at — >■ cxD. This situation changes at the cases 
of SU(4) and SU(6) at half-filling. The entropies per par- 
ticle at intermediate and strong interaction strengths are 
larger than that of the free case which is independent of 
2N. This means that for the cases of SU(4) and SU(6) 
with a larger number of fermion components, the tem- 
perature scale comparable to J is already sufficient to 
suppress spin correlations. For U/t — co, J is suppressed 
to zero, and thus T/J is set to infinity, Ssu(2N) can be 
calculated analytically as 



S. 



su{2N) _ (2A^)! 



N (A^!)2' 



(3) 



which grows monotonically as increasing 2N and satu- 
rates at In 4 as 27V — oo. 

The above picture gives rise to the main point of this 
work: the Pomeranchuk effect in SU(27V) Fermi-Hubbard 
model becomes observable at intermediate interaction 
strengths starting from the cases of SU(4) and SU(6). 
This is because the magnetic correlations are weakened 
by the large- 7V effect. We have confirmed this picture 
by performing DQMC simulations for various parame- 
ter values of T and U for the SU(6) case, which is of 
direct relevancy to the current experimental interests 
about ultracold ^^^^Yb atomsQ. The isoentropy curves 
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FIG. 2: The isoentropy curves for the half-filled SU(6) Hub- 
bard model in a 10x10 square lattice. The dashed line de- 
notes the spin superexchange scale in strong coupling regime 
J/t = 4i/U 



are displayed in Figl2l For low values of the entropy, 
say, Ssu(6)/kB < 0.7, adiabatic increasing U leads to 
a significant cooling to a temperature below the mag- 
netic superexchange scale J, which means that in this 
case, the Pomeranchuk cooling can overwhelm the de- 
creasing of J when we adiabatically increasing U, thus 
finally cool the system to the spin su per exchange tem- 
Different from Ref. 



perature regime. Different from Ref. [27| , we fix filling 
factor instead of total particle number, thus the entropy 
in the non-interacting case does not scale as TVs as they 
presented, but also scales linearly with N as the same as 
those in the half-filled Mott-insulating states. 

Next we study the simulations of charge fiuctuations 
for the half-filled S\J{2N) Hubbard model. The normal- 
ized on-site charge fluctuations are defined as 



3sti(2Af) 



TV 



(4) 



where (n^) = TV. At T oo, Ssu{2N) can be calculated 
exactly which is independent of 2TV as Ssu{2N) {T — > oo) = 
^ w 0.71, which is the upper limit. Similarly, at [/ = 0, 

5su{2N) is the same as ^ and is independent on both 
2TV and T. For the general case, we plot the DQMC 
simulation results of (5 at a relatively weak interaction 
U/t = 4 for a large range of temperatures in Fig. [3] (a). 
For all the cases, Ssu{2N) is suppressed by U from the 
upper limit of due to interaction. For the cases of 
SU(4) and SU(6), Ssu(2N){T) first falls as T increases, 
which is a reminiscence of the Pomeranchuk effect. After 
reaching a minimum at T comparable to t, dsu{2N) grows 
as increasing T. This indicates that fermions are the most 
localized at an intermediate temperature scale at which 
the spin channel contribution to entropy dominates. For 




FIG. 3: a) The on-site charge fluctuations 5 as a function of 
temperature T with parameters U/t = 4 and different values 
of 2N for a 10 X 10 lattice. The inset shows the convergence 
of 5{T) with L = 8, 10, and 12 for the SU(6) case, b) The 
normalized compressibility Hsu(2N) / i^N) v.s. T at U/t = 4. 



the SU(2) case, the non-monotonic behavior of Ssu{2N) 
is not obvious. This agrees with the picture that large 
values of 2TV enhances spin fluctuations, and enhances 
the Pomeranchuk effect. 

The physical observable compressibility k can be ex- 
pressed in terms of the global charge fluctuations as 



Ksu{2N) 



1 



1 dNf 



(5) 



where TV/ = hi is the total fermion number op- 
erator in the lattice; ^ is the chemical potential. In 
Fig. 0] (b), we plot the simulated results for the nor- 
malized Ks„(27v)/TV, i.e., the contribution to Ksu{2N) per 
fermion component. They behave similarly with each 
other. Ksu{2N) scales as 1/T like ideal gas at high tem- 
peratures, while they are suppressed at low temperatures. 
At zero temperature, Ks„(2Ar) is suppressed to zero due 
to the charge gap in the Mott-insulating states. Ksu(2N) 
reaches the maximum at an intermediate temperature 
scale which can be attributed to the charge fluctuation 
energy scale. 

Next we study spin correlations of the SU(2TV) Hub- 
bard model. The SU(2TV) generators can represented 
through fermion operators Ci^aict = 1 ~ 2TV) as Map^i = 
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FIG. 4: The normalized SU(27V) susceptibilities Xsu(2N) v.s. 
T with fixed U/t = 4 for 2N = 2, 4, 6 




FIG. 5; The normalized spin structure factor S{q) for half- 
filled SU(2iV) Hubbard model for 2N equal to (a) 2, (b) 4, 
and (c) 6. Parameter values are T/t = 0.1, and U/t — 8. 

•^I i'^l3,i ~ jfj^aisn-i- There are only {2NY ~ 1 indepen- 
dent operators due to the constraint Maa = 0. They 
satisfy the commutation relations as [Ma^^i, M~fS.j\ = 
Si.jiMas,iS'yi3 — Mjp^iSas)- Due to the Fierz identity, 
the on-site SU(2iV) fluctuations are related to the onsite 
particle number fii = J^a as J^a/i ^aii,iMi3a,i = 

(2jv+i)jv _2^^^_^)2 YVe define the SU(2iV) version 
of two point equal time spin-spin correlation as 

Mspu,{i,j) = (2jv)2 „ 1 Map,iMpa,j)- (6) 

At finite temperatures, no magnetic long-range-order 
should exist in the 2D half-filled SU(27V) model due to 
its continuous symmetry of SU(2iV). The normalized 
uniform SU(2A^) spin susceptibility is defined as 

Xsu{2N) 



The DQMC simulation results are presented in Fig. |4] 
for U/t = 4. At high temperatures, Xsu{2N) exhibits 
the standard Curie- Weiss law which scales proportional 
to 1/T. Xsu{2N) reaches the maximum at an intermedi- 
ate temperature at the scale of J below which Xsu{2N) 
is suppressed by the AF exchange. The nature of the 
ground states of half-filled SU{2N) Hubbard model re- 
mains an open question in literatures when 2A'^ is small 
but larger than 2. Nevertheless, we expect that they ei- 
ther AF long-range-ordered like the case of SU{2), or 
quantum paramagnetic with or without spin gap like in 
the large- iV limit. In either case, Xsu(2N) is suppressed to 
zero approaching zero temperature. The spin structure 
factors Ssu{2N) io) are calculated at half-filling and a low 
temperature, which are defined as 

Ssui2N) (9l - ^ E ^'"'"^M^P^nih j), (8) 

where r is the relative vector between sites i and j . The 
distributions of Ssu{2N){q} with 2A^ = 2,4,6 are plot- 
ted in Fig. El a) , b) and c) respectively. The sharpness 
of the peaks at q = (tt, tt) indicates the dominant AF 
correlations in all the cases. As increasing 2A^, peaks are 
broadened showing the weakening of the AF correlations. 

We have used the DQMC method to simulate the ther- 
modynamic properties of the two-dimensional SU(2A^) 
fermion Hubbard model at half-filling in the tempera- 
ture regime of direct interest of current experiments. The 
large numbers of fermion components enhance spin fluc- 
tuations, which facilitates the Pomeranchuk cooling to 
temperatures comparable to the superexchange energy 
scale. The local charge fluctuation, compressibility, spin- 
susceptibility, and spin structural factor are also investi- 
gated. 
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